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ABSTRACT 

We investigate the current sample of exoplanet spin-orbit measurements to determine whether a 
dominant planet migration channel can be identified, and at what con fidence. We use the predic- 
tions of Kozai migratio n plus tidal friction (Fabrvck v fc Trem ainc 2007) and planet-planet scattering 
(Nagasawa as our misalignment models, and we allow for a fraction of intrinsically aligned 

systems, explainable by disk migration. Bayesian model comparison demonstrates that the current 
sample of 32 spin-orbit measurements strongly favors a two-mode migration scenario combining planet- 
planet scattering and disk migration over a single-mode Kozai migration scenario. Our analysis in- 
dicates that between 34% and 76% of close-in planets (95% confidence) migrated via p lanet-planet 
scatte ring. Separately analyzing the subsamplc of 12 stars with T e g > 6250 K — which IWinn et all 
(2010) predict to be the only type of stars to maintain their primordial misalignments — we find that 
the data favor a single-mode scattering model over Kozai with 85% confidence. We also assess the 
number of additional hot star spin-orbit measurements that will likely be necessary to provide a more 
confident model selection, finding that an additional 20-30 measurements has a > 50% chance of 
resulting in a 95%-confident model selection, if the current model selection is correct. While we test 
only the predictions of particular Kozai and scattering migration models in this work, our methods 
may be used to test the predictions of any other spin-orbit misaligning mechanism. 



1. INTRODUCTION 

Exoplanets that transit their host stars provide oppor- 
tunity to study distant planetary systems in great detail. 
Most immediately, photometry during transit measures 
a planet's radius and density, but follow-up studies can 
provide much more. Secondary eclipse photometry, pho- 
tometric phase curve measurements, and transmission 
spectroscopy, for example, can reveal a planet's tempera- 
ture, albedo, atmospheric composition, and even weather 
patterns. While these tools investigate the physical char- 
acteristics of the planets themselves, transits also provide 
a valuable opportunity to measure details of planets' or- 
bital characteristics using the Rossiter-McLaughlin (RM) 
effect. 

The RM effect is an anomalous Doppler signal due to 
the shadow of a transiting planet crossing the face of a 
rotating star, and is measured by obtaining radial veloc- 
ity (RV) measurements of the host star during transit. 
As the approaching limb of the stellar surface is occulted, 
the total integrated radial velocity is redshifted, and as 
the receding limb is occulted, the integrated velocity is 
blue-shifted. Modeling the RM effect yields a measure 
of the angle between the orbital axis of the planet and 
the projected rotation axis of the host star, typically re- 
ferred to as A. While this angle is not itself physically 
meaningful, it constrains ip, the true angle between the 
planet's orbit and the star's rotation. In addition to be- 
ing a fundamental system parameter akin to semimajor 
axis or eccentricity, ip is a potential window into learning 
about planetary orbital migration, as different migration 
scenarios predict different distributions of tp. 

The first several RM measur ements that were made 
all indicated small values of A dWjnn et all [2001 120061: 
IWolf et aLlfeOOft INarita et al.ll2007fl . andlhus were con- 
sistent with small values of "0, which was not unex- 
pected since the orbits of all the planets in the So- 



lar System are aligned to within 7° of the solar spin 
axi s. The first mi s aligne d system, XO-3, wa s discovered 
by iHebrard et all (|2008l) , and confirmed by I Winn et al.l 
(2009a). Since then, many additional misaligned systems 
have been discovered (IWinn et al.l [2009b; Jo hnson et al.l 
120091 iTriaud et alJl201Ct ). This diversity of measured As 
suggests that planetary migration is more complicated 
than simple disk migration, which predicts planet orbits 
well aligned with stellar spins (unless the disk itself is 
misaligned), and even hints at multiple migration chan- 
nels. 

As the number of exoplanet systems with measured 
values of A increases, so does the desire and ability 
to draw conclusions based on the data. For exam- 
ple, [KbryckxAWin^ (|2009D (FW09) suggest that there 
might be two distinct populations of close-in planets — 
intrinsically aligned and intrinsically misaligned — with a 
95% probability of > 64% of planets belonging to the 
aligned population. This remarkable result was based on 
the first 10 RM measurements, which included a single 
misaligned system. 

There now exists a sample of 32 published spin-orbit 
angles, which provides a valuable opportunity to read- 
dress and extend this previous st udy, particul a rly in 
light of the recent proposition by ITriaud et al.l (|2010l ) 
that current data suggest t hat all hot Jupiters migrated 
via the Kozai mechanism (iKozail Il962t IWu et al.l 120071 : 
iFabrvckv fc Tremainel I2007I K The central goal of this 
paper is to determine whether the current sample of 
spin-orbit measurements is sufficient to begin to discern 
between the predictions of different exoplanet migration 
theories su ch as the Kozai mecha nism and planet-planet 
scattering (Nagasa wa et al.ll2008D . and if not, then how 
many more RM observations will be needed to draw more 
meaningful conclusions about the intrinsic ip distribution 
of transiting exoplanets. 

We describe our models in $2] and our analysis in S}3] 
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and <g] In $5] we repeat our a nalysis restri c ted to hot 
stars (following the suggestion of lWinn et~aT1(l20Tr3) '). In 
$6l we look to the future and ask how many RM obser- 
vations will be needed to draw more confident statistical 
conclusions. We conclude our discussion in SJ7J 

2. THE MODELS 

We test two hypotheses against each other in this pa- 
per. The first is that close-in exoplanets migrated to 
their present-day orbital locations through a combina- 
tion of the effects of Kozai cycles and tid al friction, as de- 
scribed by IHEli^x&^imimi ((2003) (FT07). Kozai 
cycles are oscillations in eccentricity and inclination of 
a close binary system caused by the presence of distant 
third companion in an inclined orbit. If a giant planet 
(assumed to have formed beyond the ice line in its pro- 
toplanetary disk) undergoes Kozai cycles that cause it 
to pass within a few stellar radii of its host star, then 
tidal friction can quickly circularize its orbit, freezing in 
a potentially large orbital inclination (ip) to the newly- 
migrated hot Jupiter. FT07 uses 1000 simulations of 
such systems, using Jupiter-mass planets with initial or- 
bital separations of 5 AU and outer 1 Mq companions 
on a 500 AU orbits highly inclined with respect to the 
planet's initial orbit. The final inclination of the resulting 
close-in planets in these simulations is their prediction of 
the hot Jupiter tp distribution resulting from this migra- 
tion mechanism (FT07, Figure 10). In the rest of the 
paper, we refer to this model as "KCTF." 

The second hypothesis is that planet-planet scattering 
is the dominan t mechanism for f ormin g close-in planets, 
as modeled by iNagasawa et all (|2008l) (N08). In their 
simulations, they study the evolution of systems of three 
Jupiter-mass planets with initial orbital separations of 
5.00, 7.25 and 9.50 AU and initial inclinations of 0.5°, 
1.0° and 1.5°. In addition to planet-planet scattering, 
they also include the effects of the Kozai mechanism and 
tidal friction, and produce a final population of close-in 
planets with a large range of orbital inclinations (N08, 
Figure 11(c)). We refer to this model as "PSTF." 

There have been many different simulati ons simi- 
lar in nature to these that we test (e.g. IWu et al.l 
(120071): IChatteriee et alJ (12001 : iFord fc Rasiol (120081) : 
Uuric fc Tremaind (j2008l) h but these are the two that 
produce the broadest distribution of ip values, and thus 
seem most likely to be able to explain the o b serve d ret- 
rograde systems. In addition, iTriaud et al.l (|2010l ) also 
compared the observed data with these two models, con- 
cluding that the KCTF model of FT07 describes the cur- 
rent data better than PSTF of N08. 

We emphasize that while the true spin-orbit angle ip 
is the physically meaningful angle, only the projected 
version of this angle A is measurable through the RM 
effect, because the inclination of the stellar rotation axis 
is unknown. If the star is rotating completely edge-on 
(/* = 90°), then A = ip, but in general the star's rota- 
tion axis may be tilted along the line of sight, resulting 
in A ip (see FW09, Fig. 3). In other words, observa- 
tion of a large value of A is firmly indicative of a large 
ip, but observation of a small A does not rule out the 
possibility of large ip. In rare cases, this unknown stel- 
lar inclination can be constrained by combining photo- 
metric determination of a rotation period with an esti- 
mate of the stellar radius and rotational line broadening 
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Fig. 1. — Probability distributions for the true (tp; top panel) 
and projected (A; bottom panel) spin-orbit angles that our two 
misalignment mechanisms produce. The KCTF tp distribution is 
taken from the simulations of final inclinatio ns of planets migrating 
through Kozai cycles + tidal damping from Fabrycky & Trcmainc 
(2007), and the planet-pl anet scattering (PSTF) prediction is taken 
from the simulations of Nagasawa et al. (2008). We use Monte 
Carlo simulations to project each i/> distribution to create the A 
distributions. Note that non-zero but small A values are preferred 
by the KCTF model, while large values (A > 60°) are preferred by 
the PSTF model. The bottom panel shows the distribution of the 
32 measured A values. 



(|Winn et al.l 120071) . Additionally, the likelihood for any 
particular transiting planet to be misaligned may be es- 
timated before even any RM measurement occurs, by 
comparing the obse rved line broaden ing to theoretical 
rotation predictions (jSchlauf man 201(J. In general how- 
ever, since the inclinations of individual stellar rotations 
are unknown, statistical methods must be employed to 
draw conclusions about ip from an observed population 
of A. 

One strategy to do this, employed by ITriaud et alJ 
(|2010f K is to statistically deproject each A measure- 
ment to form a posterior probability distribution func- 
tion (PDF) for ip for each system, assuming an isotropic 
distribution for the inclination of stellar spins, and sum 
them to form an overall ip PDF. However, this method 
has the disadvantage that isotropy of stellar spins is actu- 
ally a strong assumption, since the observed distribution 
of spin inclinations will depend on the true ip distribu- 
tion, which we are trying to determine. This is analogous 
to how the orbit inclinations of RV-detected planets may 
not always be assumed to be isot ropically distributed, as 
discussed in lHo fc Turner! (|2010D . 

Our analysis relies instead on comparing the observed 
A data directly to theoretical predictions of the models, 
in A-space. This requires that we determine a probabil- 
ity density function (PDF) for A corresponding to each 
of the ip distributions we wish to test. Since both of our 
ip predictions are the results of complicated simulations, 
we use a Monte Carlo procedure to perform this trans- 
formation. 

We fix an observer-oriented spherical coordinate sys- 
tem and simulate a large number of systems as fol- 
lows. First, we populate 10 6 stars with transiting plan- 
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ets with orbital inclination vectors O assigned according 
to the distribution of known transiting planet inclina- 
tions. Then the stellar spins S are assigned relative to 
the planet population, according to the predicted ip dis- 
tributions of KCTF and PSTF. This is accomplished by 
selecting a ip from the distribution implied by the simu- 
lations of FT07 or N08, treating this if> as a polar angle 
relative to O, and assigning S to have an azimuthal an- 
gle around O, chosen uniformly on (0, 2n). Then S is 
transformed back into the observer-oriented coordinate 
system, in which A is simply the azimuthal angle of S 
about the line of sight. The probability distribution for 
A is then determined from the distribution of these re- 
sulting A values. Figure [T] illustrates the original ip and 
derived A distributions for both misalignment models, as 
well as the current A data. 

In addition to comparing these two misalignment mod- 
els, we also consider the possibility that there might be 
multiple migration mechanisms, inspired by the conclu- 
sions of FW09, who found the data available at the time 
(10 systems) favored a model with planets drawn from 
two distinct distributions: one perfectly aligned (tp = 0) 
and one isotropically misaligned. After introducing our 
methods in |£3j we explore in |J4] what we can learn if 
we assume only a fraction / of systems are misaligned 
according to one of the above mechanisms, with the re- 
maining 1 — / fraction being perfectly aligned. 

3. DATA ANALYSIS 

The goal of our analysis is to select the model that 
describes the data best, and to determine the confidence 
with which we can make that selection. We do this first 
assuming that all planets are misaligned, and then again 
allowing for an aligned population. The following sub- 
sections outline the details of these steps. The data we 
us e are the 28 RM measurements compiled in Table 1 
of|W inn et all (|2010f ) (with 5 of these systems updated; 
Winn, in prep), plus HAT-P-14 (185° ± 4.5), (Winn, in 
prep), HAT-P-4 (- 15° ± 16) (Winn, in prep), and XO-4 
(jNarita et afllMtCl) . A summary of the results of all our 
analysis is in Table 1. 

3.1. Which misalignment mechanism is preferred? 

We use the Bayes factor (e.g. ?) for our model selec- 
tion statistic, which in the simple case of comparing two 
models with no free parameters reduces to a likelihood 
ratio, as follows: 



n = io gl 



£koz({A}) 

^scat ({A}) 



(1) 



where £koz({A}) and £ SC at({A}) are likelihoods of the 
observed data {A} under the two different misalignment 
models. K > favors the KCTF model, and TZ < 
favors the PSTF model. The likelihoods are calculated 
as follows: 

N 

C M ({\}) = l[PM(Xi), (2) 



where 



-Px(Ai) 



ISO 



(3) 



with Pi(X) being the probability distribution of the ith A 
measurement (which we take to be a Gaussian centered 
at the published value with width as the published error 
bar), and pm,\W being the A PDF for the model in 
question (A4 — {Koz,scat}). 

Using the current set of 32 A measurements, we calcu- 
late TZ = 1.34, which favors KCTF over the PSTF. The 
following section explains how we quantify the strength 
of this model selection. 

3.2. Confidence Assessment 

Within the dichotomous paradigm of comparing two 
misalignment models we can determine the confidence in 
our model selection by answering the following question: 
"Given a measured value of TZ = 1.34, which favors the 
KCTF model, what is the probability that the KCTF 
model is actually correct?" Or, more generally, given any 
measured value of TZ, how can we quantify the confidence 
in the implied model selection? 

To address this question, we perform the analysis de- 
scribed above on many simulated data sets. Starting 
with 10 5 randomly-generated spin-orbit systems for each 
misalignment model, constructed according to the proce- 
dure described in we randomly select 32 A values from 
each underlying model. Each simulated A measurement 
is the exact value of A drawn from the simulations, per- 
turb ed by a measurement e rror o~\, assigned using Eq. 16 
from lGaudi fc Winnl (12001 : 



1 



v sin IWN 7 



(1 - fe 2 )cos 2 A + 36 2 



6 2 (1 - 6 2 ) 



1/2 



. ( 4 ) 

where a is the single-point RV measurement uncertainty, 
usin/* is the projected rotational velocity of the star, N 
is the number of RV points in transit, 7 is the planet-star 
radius ratio, and b is the transit impact parameter. For 
these simulated data sets, we randomly assign v sin/*, 7, 
and b by drawing randomly with replacement from the 
current set of all transiting planet system^. We take 
(i = 5 m/s and TV = 20 for each simulated measurement. 

We repeat this data simulation procedure 5000 times 
and calculate TZ (Eqn. [T]) for each data set, giving us 
an understanding of the expected distribution of 1Z if ei- 
ther of these models do describe the actual underlying ip. 
Using these simulations to construct PDFs for 7Z under 
each model (pn,Koz(TZ) and Pn,sca,t{TZ)) we can then ask 
what the probability is of either model being true, given 
a measurement of 1Z. Applying Bayes' theorem with a 
uniform prior on which model should be correct, we may 
write: 



Fr(M\n) = ^ M) = , (5) 



^Pr(^|M) ^PnMK)' 



M 



M 



where again M — {Koz, scat}. For our specific case, this 
becomes: 



Pr(Koz|ft = 1.34) 



PTC,Koz(l-34) 



PK,Koz(l-34) +PTC,scat(l-34) 



= 0.96, 

(6) 



1 as listed at www.exoplanets.org 
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Fig. 2. — The relationship between our model selection statistic 
TZ, or the logarithm of the likelihood ratio, and a model selection 
confidence. We randomly draw 5000 sets of 32 A values from each 
model and measure TZ for each of these data sets. The top panel 
show the distribution of R values attained from these data simu- 
lations. The confidence in the model selection at any particular 
TZ (bottom panel) is determined by the relative heights of the two 
distributions at that TZ. The current data strongly favor the KCTF 
model. 



giving 96% confidence in the KCTF model. Figure [2] 
illustrates this confidence-assessment procedure. 

This resul t app ears to support the conclusion of 
iTriaud et al.l (|2010f) , who claim that the current RM data 
points to the FT07 KCTF model as explaining the forma- 
tion of hot Jupiters better than the N08 model. However, 
as there are reasons to be skeptical that the Kozai mech- 
anism could plausibly be responsibl e for the form ation 
of all hot J upiters, both theo retical ()Wu et al.ll2007f) and 
empirical (Schlaufman 2010), we take our analysis one 
step further and consider what we may conclude if we al- 
low for two distinct populations of systems: aligned and 
misaligned. 

4. MULTIPLE MIGRATION CHANNELS? 

Of the 32 A measurements to date, 14 have |A| < 10 
and 19 are within 2a of A = 0. Given that KCTF pre- 
dicts many more systems with small A than does the 
PSTF model, it is thus not surprising that it is pre- 
ferred over PSTF in the above analysis. But what if 
only a fraction of planets migrate via a misaligning mech- 
anism while th e rest migrate t hrough a process such as 
disk migration (|Lin et al.lll996l) that preserves spin-orbit 
alignment? Analyzing the first 10 A measurements (that 
included only a single significantly misaligned system), 
FW09 found such a two-population model to be their 
best selection. Especially given the difficulties of ex- 
plaining all hot Jupiter mig ration usi ng the Kozai mecha - 
nism alone (|Wu et al.ll2007l) . and that lSchlaufmanl (|2010f) 



found that fewer systems seem to misaligned than would 
be predicted based on this prediction, it seems prudent 
to investigate how allowing for an aligned population af- 
fects conclusions about the intrinsic tp distribution. 

Accordingly, we add a parameter to each of our models: 
/, the fraction of systems that are misaligned, with ifj — 
(and thus A = 0). This requires us to modify our model 
selection procedure Q]4.2[) . but it also enables us to ask 
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Fig. 3. — The posterior distribution of / for each of our two- 
migration-channel models, conditioned on 32 observed A values. 
We show that if the true misaligned population were misaligned 
according to the KCTF distribution, then we expect that about 
64% of close-in planets are misaligned. Similarly, if the misaligning 
mechanism were PSTF, then the most likely value of / is close 
to 55%. Strictly speaking, this figure shows the likelihood of the 
current data as a function of / for each model — if these curves 
were normalized they would be proper probability distributions, 
given our flat prior for /. However, plotting them in un-normalizcd 
form is illustrative, since we use the ratio of the areas under these 
curves (the marginal likelihood) as our two-mode model selection 
criterion. 



an additional question: 

/ (SU? 



what do the models imply about 



4.1. What fraction of systems are misaligned? 

We may use Bayes' theorem to calculate the probabil- 
ity distribution for the misaligned fraction / for each of 
our two models, conditioned on the 32 observed A values. 
For this particular case, we may write Bayes' theorem as 
follows: 



PfMf\W) = 



C M ({X}\f)p(f) 
C M ({X}\f)p(f)df 



(7) 



where as before {A} is the set of observed A values and 
M. represents a particular model (either 'Koz' or 'scat'), 
P/,.m(/|{A}) is the posterior probability distribution for 
/ under model A4, £wi({A}|/) is the likelihood of the 
data given a particular / under M, and p(f) is the prior 
probability distribution for /, which we take to be flat 
between and 1. The denominator is the normalization 
factor, also known as the marginal likelihood. The poste- 
rior probability distribution for / allows us to infer con- 
clusions about / for a particular model given the current 
data. The likelihood function is calculated the same way 
as Eq. [21 except that now the probability distribution for 
A is dependent on /: 



Px, M (Mf) = fxpMW + (i-f)s(x), 



(8) 



where pmW are the A distributions we calculated in 
(bottom panel of Figure [l}, and 5(A) is the Dirac delta 
function. 

Fig. [3] illustrates £x({A}|/) as a function of / for each 
model. We measure the most likely values and their 
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Fig. 4. — The relationship between our model selection statistic 
and model selection confidence, for comparing our two-mode mi- 
gration models, where only a fraction / of planets are misaligned 
and the rest are in perfectly-aligned systems. 1Z in this case is 
the logarithm of the ratio of the marginal likelihoods of the two 
models, and the 1Z distributions are generated by measuring the 
1Z values for data simulations at different /. The PSTF + align ed 
model is preferred over the KCTF + aligned model. See i|4,2l for 
details. 

symmetric 95% confidence ranges for / by normalizing 
the likelihoods and computing the cumulative distrubu- 
tion functions. If the KCTF + aligned model is cor- 
rect, / lies between 0.40 and 0.87 with 95% confidence, 
with the most likely value being 0.64. Similarly, if the 
PSTF + aligned model is correct, / lies between 0.34 
and 0.76 with 95% confidence, with the most likely value 
being 0.55. Both models indicate a significant fraction 
of aligned systems, with the PSTF model indicating that 
nearly half of systems might be intrinsically aligned. So 
the next question to ask is in this two-channel picture of 
planet migration, which model do the data favor? 

4.2. Two-mode model selection and confidence 
assessment 

Since the models we are now comparing each have an 
unknown parameter /, we redefine our model comparison 
statistic to be the ratio of the marginal likelihoods (the 
denominator of Eqn. or the area under the curves in 
Figure [3|) of the two models: 



/ f 1 



1Z = loa 



V7o 



£ Koz ({\}\f) P (f)df 



({\}\f) P (f)df 



(9) 



where again 1Z > favors KCTF + aligned and TZ < 
favors PSTF + aligned. This time, we calculate TZ = 
-1.16, which favors PSTF over KCTF. 

Confidence assessment requires data simulations as be- 
fore, except now we simulate data sets for each model on 
a grid of / values, to determine the behavior of TZ as a 
function both of misalignment model and /. So for each 
of 20 equally-spaced values of / between and 1 , we ran- 
domly draw 1000 sets of 30 systems, as described in ^3-21 
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Fig. 5. — The relationship between our model selection statistic 
and model selection confidence, comparing the preferred single- 
mode migration scenario (KCTF) to the preferred two-mode model 
(PSTF + aligned). PSTF + aligned is strongly preferred over 
single-mode KCTF. 



simulating the aligned fraction of planets by giving each 
simulated A a probability (1 — /) to be re-assigned to 
A = before being perturbed by the measurement error. 

These simulated 1Z values give us an 1Z PDF for each 
/ for each model — essentially empirical two-dimensional 
likelihood functions: £tc,Koz(7^, /) and C 

(ft,/), 

where 1Z is continuous but / is only sampled at 20 
points between and 1. Since we already calculated the 
posterior probability distribution of / given the current 
data for each of our models above, we may marginalize 
these likelihood functions over / to calculate a properly 
weighted one-dimensional PDF for 1Z: 

20 

PnMK) =X^,x(ft/ i )Pr(/ i |{A},>f). (10) 

i=l 

The term Pr(/i|{A}, Ai) is the probability obtained 
by integrating p/ i: ^((/|{A}) (Eqn. [TJ Figure|3J) between 
/, — A//2 and fi + A//2, where A/ is the spacing be- 
tween successive / values in our simulations. Using these 
likelihood functions, we may calculate the probability of 
either of our models being true, analogously to Eqns [5] 
and[H We find that our measurement of 1Z = —1.16 gives 
a confidence of 94% for the PSTF model, illustrated in 
Figure |3J 

Thus, while the current sample of A measurements ap- 
pears at first to be evidence for Kozai migration as pre- 
dicted by FT07, the model preference changes in favor 
of the N08 scattering model if allowance is made for the 
existence of some systems forming hot Jupiters via some 
mechanism that preserves alignment. And since what 
current RM measurements can tell us about planet mi- 
gration mechanisms changes significantly depending on 
whether we assume one or two migration channels, we 
next explore whether current data allows us to distin- 
guish between those scenarios. 

4.3. One channel or two channels? 

We approach this question the same way we have hith- 
erto approached the other model selection questions: de- 
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Fig. 6. — The relationship between our model selection statistic 
and model selection confidence, but using only the 12 hot stars 
(jfcff > 6250) as our data. One-channel migration favors the PSTF 
mechanism for this subset of the data, but not conclusively. The 
1Z thresholds required to reach 95% confidence in either model for 
this sample size are marked. 
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Fig. 7. — The posterior distribution of / for each of our two- 
migration-channel models, conditioned on the 12 observed hot star 
(T cff > 6250) A values. As lWinn et~aLl H2010I ) point out, hot stars 
appear to be misaligned more often than the overall stellar popu- 
lation. Because / = 1 has a high probability under both models, 
we do not further pursue two-mode migration model comparison 
for this subsample. 



fine a model selection statistic TZ, calculate TZ for the 
current data, and determine confidence based on data 
simulations. Our model selection statistic in this case 
becomes: 

n = io gl0 , (ii) 

where £ onc is one-mode likelihood as used in Eq. [TJ and 
£two is a marginalized likelihood as used in Eq. Since 
KCTF is preferred for one-mode migration and PSTF 
is the preferred two-mode model, we take C onc to be 
the likelihood of the data under the KCTF model and 
£two to be the marginal likelihood of the data under the 
PSTF+aligned model. We calculate TZ = -3.00, which 
strongly favors the PSTF + aligned model over the one- 
mode KCTF model, with a confidence of 99% (Figure 

5. HOT STARS 

We have shown that the KCTF prediction for the dis- 
tribution of ip for migrated planets does not adequately 
explain the observed distribution of A, and that current 
data favor a combination of well aligned systems and 
systems misaligned via planet-planet scattering. How- 
ever, it may be that the current population we observe 
is not representative of the initial misalignm ent distri- 
bution. For example, iMatsumura et all ([2010) conclude 
that tidal effects may be import ant in damping out initial 
misalignment in some systems. iWinn et al.l (|2010| ) note 
that most of the misaligned planets that have been ob- 
served to date are around stars hotter than 6250 K. Based 
on this empirical finding, they suggest that perhaps all 
HJs migrate through some misaligning mechanism, and 
that cooler stars with deep convective zones have their 
envelopes tidally torqued into alignment by the planet on 
a relatively quick timescale, thereby erasing the evidence 
of the initial misalignment. If this were indeed the case, 
then an important key to understanding HJ migration 
would be spin-orbit measurements of planets transiting 
hot stars, since these systems would presumably have 



maintained their primordial misalignment. 

Following this line of inquiry, we repeat the analyses of 
§3 and §4 restricted to only the current sample of 12 hot 
stars (T c ff > 6250 K). Of our two one- mode migration 
models, PSTF is favored over KCTF, with 82% confi- 
dence (TZ = —0.73; Figure [6]). We also explore the two- 
mode migration hypothesis, allowing for hot stars to also 
have an intrinsically aligned population. The posterior 
distributions for / under the two models conditioned on 
the hot star data are illustrated in Figure [71 showing 
that both models favor almost complete misalignment. 
In fact, one-mode vs. two-mode model selection for this 
subset of the data gives TZ = —0.33, favoring the single- 
mode hypothesis. For this reason we do not pursue any 
further the idea that there exists an intrinsically aligned 
population among hot stars, as the current data do not 
merit the additional model complication. 

Thus, if the distribution of hot star A values represents 
the primordial alignment distribution for all stars, the 
current 12 A observations hint that the PSTF prediction 
of N08 describes planet migration better than the KCTF 
model of FT07. However, there are not yet enough data 
for this model selection to be conclusive. In the next 
section, we discuss how many more RM measurements 
will likely be needed in order to make hot star model 
selections confident. 

6. HOW MANY HOT STAR RM OBSERVATIONS 
ARE NEEDED? 

With the full sample of 32 RM observations to date, 
we are able to confidently state that a combination of 
well aligned systems and systems misaligned via planet- 
planet scattering (the PSTF model) explains the current 
data better than migration via the Kozai mechanism and 
tidal friction alone (KCTF). Only considering the sam- 
ple of 12 hot stars, our present conclusions are w eaker. 
Inspired by the work of iSwift fc Beaumo nt (2010), who 
calculated the sample size of protostellar cores necessary 
to reliably distinguish a power law from a lognormal dis- 
tribution, we wish to quantify how many additional hot 
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Total number of hot star systems with measured A 

Fig. 8. — The probability with which a given-sized future sample 
of hot-star (T e g > 6250 K) A values (in addition to the 12 that 
currently exist) will result in a confident selection (> 95%) of one 
of our misalignment models over another. We use data simulations 
of various sample sizes (described in [|6j to define the confidence 
thresholds and to determine how likely those thresholds are to be 
reached for each particular sample size. For example, with 28 more 
measurements (40 total), we have a 75% chance of selecting the 
PSTF model at > 95% confidence. 

star A measurements will likely be needed in order to 
measure 1Z values indicative of confident (> 95%) model 
selection between the KCTF and PSTF models. 

This requires determining two things: first what values 
of 7Z correspond to 95% confident model selections for a 
given sample size N, and secondly how likely we are to 
actually measure such a confident value once we have N 
observations in hand, given the current set of 12 as our 
starting point. 

Both questions can be addressed by data simulations. 
The TZ thresholds that represent 95% confidence for a 
given N can be determined by simulating many data 
sets of sample size N and using Eq. [5] to calculate model 
selection confidence as a function of 1Z. This has already 
been illustrated for N = 32 in Figures [21 01 and [5] and for 
N = 12 in Figure [5] We repeat this procedure for dif- 
ferent sample sizes, defining 1Z thresholds for our model 
selection for a series of N values up to N — 100. 

We also use data simulations to determine the proba- 
bility of a future experiment resulting in a confident 1Z 
measurement. To do this we again simulate many data 
sets of the same sizes as above and measure the 1Z dis- 
tributions of the experiments, but this time the first 12 
measurements in each data set are fixed to be the current 
hot star measurements. Then we use the 1Z distributions 
to predict how often threshold values are reached for each 
model. This probability is plotted against sample size in 
Figure [5] 

We learn from this that in order to have a good chance 
of obtaining a hot star data set that confidently selects 
either of our single-mode migration models, we will likely 
need a total data set of about 80-100 hot star RM obser- 
vations (Figure H]). More optimistically, there is > 50% 
chance of confident model selection with only a factor 
of ^3 more observations (40 total) if the PSTF model 
is correct. Given the pace at which this field is grow- 
ing, this may be reasonably expected to happen within 



the next few years. On the other hand, if the KCTF 
model (or some other model) is a better description of 
reality, then it will likely take more observations to de- 
termine this, given the preference of the current data for 
the PSTF model. 

7. DISCUSSION AND CONCLUSIONS 

The study of exoplanet spin-orbit angles is advancing 
rapidly, with 32 measured projected spin-orbit angles and 
more to come as ground- and space-based surveys con- 
tinue to detect transiting planets. We have analyzed the 
current sample and quantify what may be inferred about 
the distribution of true spin-orbit angles ip. In particular, 
we ask whether the current data are sufficient to test the 
predicted distribution of ip from specific migration mech- 
anisms, u sing the Kozai cycles + tidal friction (KCTF) 
model of iFabrvckv fc Tremaind (2007) and the p lanet- 
planet scattering model of iNagasawa et aD (|2008l ) (that 
also includes the Kozai effect and tidal friction; PSTF) 
as test models. 

We find that conclusions about which migration mech- 
anism is responsible for misalignment depend crucially 
on the assumption of whether there exists a popula- 
tion of intrinsically aligned systems (ip — 0). Without 
allowing for an intrinsically aligned population we find 
that the KCTF model is favored over the PSTF mecha- 
nism ($3]), but allowing for this population we find that 
PSTF is favored ($IJ, with the most likely fraction of 
misaligned systems being 0.55, and between 0.34 and 
0.76 with 95% confidence. We also find that this two- 
mode migration model (PSTF + aligned) is significantly 
favor ed over sing l e-mod e KCTF migration. This agrees 
with iSchlaufmani (|2010D . who also concluded that there 
is likely to be an aligned population, based on detecting 
fewer likely-misaligned systems than expected based on 
predictions of the KCTF model. 

These results may be an indication of two migration 
channels for close-in planets, one that acts gently, pre- 
serving the alignment of planet orbits with the stellar 
spin, and one that acts impulsively, causing misalign- 
ment. The gentle mec hanism might well be disk mi- 
gration (|Lin et al.lll996| ). and our analysis suggests that 
the misaligning mechanism is not solely the Kozai ef- 
fect but rather some mechanism that distributes ifi more 
broadly, such as planet-planet scattering in combination 
with the Kozai effect. Th is accords with the conclusion of 
iMatsumura et all (|2010T ) that some of the close-in plan- 
ets with non-zero orbital eccentricities are likely to have 
been formed by planet-planet scattering and subsequent 
tidal circularization. 

Focusing on the subsampl e of 12 hot star (T B ff > 6250 
K) A measurements, which iWinn et al.l (|2010f ) predict 
might be the only systems to maintain their primordial 
misalignments, we find that the data prefer the single- 
mode PSTF model over KCTF, with a confidence of 85%. 
We also find that a single migration mechanism is suf- 
ficient to describe the current hot star A distribution, 
without including an intrinsically well aligned fraction. 

Looking to the future, we also calculate the number of 
additional hot star A measurements necessary to achieve 
> 95% confidence in the hot star model selection (SjB|). 
We find that if either of our single-mode mechanisms 
describes the tp distribution around hot stars, then a 
total data set of about 80-100 A measurements should 



TABLE 1 

Results summary 



Data 


Model 1 


Model 2 


/ a (model 1) 


/ a (model 2) 


n h 


Confidence 


All 32 systems 


KCTF d 


PSTF C 






1.34 


96% 


All 32 systems 


KCTF + aligned 


PSTF + aligned 


0.64 ±0.24 


0.55 ±0.21 


-1.16 


93% 


All 32 systems 


KCTF 


PSTF + aligned 




0.55 ±0.21 


-3.00 


99% 


12 hot systems' 


KCTF 


PSTF 






-0.73 


81% 



a maximum-likelihood value, with symmetric 95% confidence range; only applicable to two-mode models 
b 1Z > favors model 1; 1Z < favors model 2 

c Degree of belief that the model selection is correct; based on Monte Carlo simulations (§3.2) 

d Kozai cycles + tidal friction (Fabrvckv & Trcmainc 2003) 

c Planet-planet scattering + Kozai cycles and tidal friction (Naeasawa ct al. 2008) 
f host star T cff > 6250 K 



definitely be sufficient to solidify which is the preferred 
model, with a > 50% probability of confident model se- 
lection with a total data set of only about 40, if scattering 
is indeed the best explanation for the ip distribution of 
close-in planets. 

Thus, we suggest that if RM studies wish to be able 
to identify migration mechanisms through measuring A 
distributions, they should concentrate on measuring A 
for planets around hot stars. Of course it is conceivable 
that migration mechanisms themselves might be different 
around different types of stars, in which case hot star A 
measurements might not tell us anything about cool star 
migration, but that is the kind of question that will only 
be able to be explored when much more data is available. 

Finally, we emphasize that the results in this paper 
are illustratory more than definitive, as we have only 
tested only two very specific misalignment models. Con- 
sequently, we encourage continued theoretical work pre- 
dicting ip distributions, as the analysis we present may 
be applied to any prediction. We can use the results of 
this paper as a guide for what to expect from such future 
analyses. For example, models that favor larger values of 
A more strongly than the present KCTF prediction (as 
does the PSTF model) are likely to be preferred. We also 
suggest that an interesting question to pursue would be 
self-consistent planet formation and dynamical evolution 
calculations, in order to explore planet-planet scattering 
in the context of realistic formation sc enarios (in con- 
trast t o the fixed initial conditions of thelNagasaw a et al.l 
(2008) simulations. A particularly intriguing angle to ex- 
plore in such work would be whether a plausible expla- 
nation for the observed trend in A with stellar temper- 
ature/mass might be explainable by more massive stars 



tending to form more massive planets in closer proximity, 
so as to make planet-planet scattering (and thus spin- 
orbit misalignment) more common among earlier-type 
stars. Migration models might also be combined with 
models that predict that large values of ip might origi- 
nate from the host s tar itself being tilted relative to the 
disk (|Lai et al.H2010f ). rather than solely from the effects 
of planet migration, though current observations suggest 
that p rotoplanetary disks tend to be aligned with stellar 
spins (| Watson et al.ll2010h . 

Planet migration has been a mystery ever since the 
first hot Jupiters were discovered, with very little ob- 
servational evidence to substantiate theories. As fossil 
remnants of planet migration, spin-orbit angles are key 
to understanding the origins of these close-in planets, 
and the first few years' worth of A measurements are be- 
ginning to give substantial clues. Many more transiting 
planets will be discovered in the near future thanks to the 
productivity of transit surveys, and as long as RM mea- 
surements of these systems continue, our understanding 
of planet migration will continue to improve. 
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ture and content of this paper from Josh Winn and Dan 
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course of this project, and also acknowledges the Penn 
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ify his thinking about model selection. J. A. J. thanks 
Jon Swift, Michael Cushing, Brendan Bowler, Justin 
Crepp, and Ed Turner for illuminating discussions over 
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methods. 
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